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Abstract 

A scheme based on a unifying q-deformed algebra and associated with a generahzed Lax operator 
is proposed for generating integrable quantum and statistical models. As important applications 
we derive known as well as novel quantum models and obtain new series of vertex models related 
to q-spin, q-boson and their hybrid combinations. Generic q, q roots of unity and q ~f I yield 
different classes of integrable models. Exact solutions through algebraic Bethe ansatz is formulated 
for all models in a unified way. 
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1 Introduction 

Integrable quantum systems in 1 + 1-dimensions and statistical models in 2-dimensions belong to an 
exclusive and important class of models giving exact results. However there seems to be no established 
scheme for generating all such models along with their Lax operators and i?-matrices in a unified way. 
Our aim therefore is an effort in that direction for constructing a wide range of known as well as new 
classes of models putting particular emphasis on statistical models. 

The integrable systems may be defined by the property that they possess N-number of conserved 
quantities, which are independent and commute among themselves: [cn,Cm] = 0, where N is to 
the degree of freedom of the system. Such integrability also leads to the exact solvability of the 
model, since taking conserved quantities as action variables one can adopt the action-angle variable 
description. For managing conveniently such a rich structure we introduce first a generating function 
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r(A) depending on some extra parameter A, known as the spectral parameter, such that one can recover 
the infinite number of conserved quantities including the Hamiltonian as the expansion coefficients 
of Inr(A) = ^jCjX^. The crucial integrability condition may then be defined in a compact form 
as [r(A),r(/i)] = 0, from which the commutativity of Cj's follows immediately by comparing the 
coefficients of different powers of A, /i. The partition function Z of the related statistical model on the 
other hand can be constructed from t(A) as Z = tr(r(A)^). 

To be convinced of the amazingly wide variety and range of the integrable models, let us name a 
few of them [||, which allow both classical and quantum variants. They are field and lattice models 
of relativistic or anisotropic type such as sine-Gordon field and exact lattice |2[ model and similarly 
Liouville model, quantum derivative nonlinear Schrodinger (DNLS) equation Q together with their 
exact lattice versions, quantum spin model like XX Z chain relativistic quantum Toda chain Q etc. 
Models can also be of nonrelativistic or isotropic type such as NLS model ^ and its lattice variant, 
XXX spin chain Q, Toda chain etc. Among the statistical integrable models the most simple 
and well known is the 6-vertex model This model is defined on a square lattice with a random 
direction on each bond ( left or right on horizontal, up or down on vertical) (FIG. la)) constrained 
by the ice rule, that the number of incoming and outgoing arrows at each vertex are the same. This 
leaves only 6 possible configurations giving that many Boltzmann weights (BW) uJij-k,i at each vertex 
point, coining the name of the model. These 6 weights may be given by the nontrivial matrix elements 
of the i?-matrix (|2.2|) and the partition function is expressed as Z = J2configI\a,b,j,k^a,j;b,k{X)- 

2 Lax operator approach 

For describing the intricate integrable structures mentioned above one can not however start just from 
the Hamiltonian of the model as conventional in physics, since now the Hamiltonian is merely one 
among many commuting conserved charges. We on the other hand have to allow certain abstractions 
and start with an unusual type of matrix called the Lax operator Lj{X) defined at each site j of a 
1-dimensional discretized lattice. The matrix elements of the Lax operator, unlike in an usual matrix, 
are themselves operators acting on some Hilbert space. For ensuring the integrability of a system 
the related Lj{X) generally should satisfy certain commutation relations given by the ultralocality 
condition Lj{X) Lk{ii) = (/ ® Lk{ii)) {Lj(X) <S> I) at different lattice points and the Yang-Baxter 
equation 

R{X - fi)Lj{X) Ljifi) = {10 Ljin)) {Lj{X) I)R{X - fi). (2.1) 

at the same points j = 1,2, . . . , N. The R{X — /i)-matrix has spectral parameter dependent c-number 
elements, which we take as a 4 x 4 matrix with 

Rll = i?22 = sma{X + 1), Rll = R^ = sina(A), R^ = RI\ = sin a. (2.2) 

Since our intension is to establish the integrability which is a global property, we have to switch 
from this local picture at each site j to some global one by defining a matrix T(A) = Y\^=iLj{X). 
Multiplying therefore the YBE ( p.lj) for j = 1,2, . . . ,N and thanks to the ultralocality condition one 
arrives at the global YBE: i?i2(A - ^) r(A) ® T{fi) = {I (g> T{^i) T®I{X) RuiX - /i) having 
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exactly the same structural form. Invariance of the same algebraic form also for the tensor product 
of the algebras as revealed here exhibits an underlying deep Hopf algebra structure. Defining further 
r(A) = trT(A), taking trace from both sides of the global YBE and canceling the i?-matrices due to 
the cyclic rotation of matrices under the trace we arrive finally at the trace identity through t(A) 
defining the quantum integrability of the system. 

Therefore we see that the integrable systems may be defined by their Lax operators satisfying 
the ultralocality condition and the YBE together with the associated i?-matrix. Note that all the 
varied integrable models listed above have their representative L operators as 2 x 2 matrices, but with 
diverse forms expressed through different basic operators like spins a, bosons ipi V'^ oi' I'sal canonical 
operators u,p etc. However surprisingly, the i?-matrix associated with all of them appears to be the 
same and given by the known trigonometric solution ( p. 21) as for the 6- vertex model (or its q = e*" — > 1 
, i.e rational limit). Therefore a natural question arises asking whether there can be a unified model 
associated with ( |2.2| ) such that all above models can be derived from it in a systematic way. 



3 Unified model 

We show that it is indeed possible to find such an integrable unified model and justify consequently the 
sharing of the same i?-matrix by all other integrable models, which will be obtained through various 
reductions and realizations of this single model. Such a general model must be a quantum, lattice 
model, which is relativistic or anisotropic and extendible to be inhomogeneous with inbuilt quantum 
parameter h, lattice constant A, deformation parameter q and the inhomogeneity parameters {c}), 
such that it can generate lattice/ field models, quantum/ classical models, relativistic (anisotropic)/ 
nonrelativistic (isotropic) variants and similarly inhomogeneous/homogeneous models at the corre- 
sponding limits of these parameters. 

Since the integrable models can be represented by their Lax operators, we propose the defining 
form for our unified model as 

L\ \0=[ ^^,^_,3^^_,^_^,3 j, <i = e,i = e ,. = 2sma. (3.3) 

with the basic operators satisfying the quadratic algebra 

[5^ S^] = ±5^, [5+, S-] = (m+ sm{2aS^) + M" cos{2aS^)) ^7^, [M^, •] = 0, (3.4) 

where are the central elements expressed as = ztV ibl(c^C2 ± c^c^) through the commuting 
set of elements c^. The algebra ( |3.4D is a novel q-deformed algebra and includes known q-spin as well 
as q-boson algebras as particular reductions. The quantum integrability of this model and hence all 
other models derived below from it is guaranteed, since ( |3.3| ) with ( |3.4| ) associated with i?-matrix (|2.2| ) 
satisfy YBE (pi)). 

It is crucial to note that we would define the Boltzmann weights (BW) of our vertex models not 
by the i?-matrix as conventional, but through the elements of the Lax operator (|3.3|) as 1?^(X) = 
'^a,j;b,A:(A) by using matrix representations of the unifying q-deformed algebra (^.41) . 



3 



We may find an important representation of tliis algebra tlirough canonical fields as 

S^ = u, S+ = e-'Pgs{u), S-=gsiu)e'P. (3.5) 

with operator function 

1 I 

gs{u) = ( K + sin a{s - u){M+ sin a{u + s + I) + M" cos a(n + s + 1))) 2 (3.6) 

sin a 

containing extra free parameters k and s. The unified model represented by ( |3.3| ) is a quantum 
integrable model and may be considered as a generalized lattice SG model through realization ( |3.5D ,( 
|3.6| ). We are now in a position to generate the whole range of integrable models, known as well as new, 
through various choices of the central elements as well as by mapping into different realizations. 
In fact ( ^^ ) can be directly mapped through the spin-i operators a or realized further in bosonic 
operators ip^'^^ as ip = e~*^((s — u))^,N = s — u, etc. We may repeat similarly the whole construction 
at limit a ^ for generating the rational cases, or consider q as roots of unity to get the restricted 
models. On the other hand at A —> we recover the field models, while h ^ yields as usual 
the corresponding classical dynamical systems. As we show below, the choice for the inhomogeneity 
parameters can be of two types: either as constant parameters reproducing generally known models 
or as site j dependent functions, which generates new classes of inhomogeneous or hybrid integrable 
models. 

Motivated by the form of ( |3.5| ) we find also a matrix representation of the unifying algebra (|3.4| ) 

as 

< S,'m\S^\m,S >= rn6rn,m, < 3,71118^1771,3 >= f^{m)6m±l,m, (3.7) 
with /^(m) = f^{m + 1) = gs{m) and using it construct the BW for our unified vertex model as 

u:±,k;±,k{\) = cie^"^^^'") + cie-^"^^^"^), u:+,k;-,k^i = c^-,fe-i;+,fc = 2g.(fc - 1) sin a, m = s + 1 - k, 

(3.8) 

where k G [1, -D], depends on the dimension D of the matrix-representation of the q-algebras. Possible 
reductions of this unified model ( |3.8[ ) would yield new series of vertex models (see FIG.l) with the 
familiar ice-rule is generalized here as the color conservation a-\-j = h + k for determining the nonzero 
BW. 

The eigenvalue solution of the transfer matrix related to the unified model can be found exactly 
through the algebraic Bethe ansatz, which therefore would give also exact solutions for all other 
quantum and statistical models we construct here in a unifying way. 

4 Construction of integrable models 

4.1 Trigonometric class with generic q 

1. The simplest constant choice = giving M~ = 0, M"*" = 2, reduces ( |3.4D to the well 
known Uq{su{2)) q-spin algebra |] [S^,S^] = ±S^,[S+,S-] = [2S% = and (|3j) to the 
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corresponding Lax operator 

\x + s\ s 



L,sM>^)={ '"^l^ j^_^3^ ), (4.9) 



expressed through the q-spin. The simplest representation S = ■^a reduces (O) further to ( |2.2[) and 

1 

"2 



recovers the weh known XXZ spin-h chain and also the 6-vertex model as the related statistical 



system. 

On the other hand, for canonical representation (|3.5| ) the form (3.6) reduces to 
g{u) = 2 Q, [1 + cos a{2u + 1)] 2 recovering the known integrable lattice sine-Gordon model. However 
for the corresponding statistical model the BW reduced from (|3.8| ) takes the form uj±j-±j{u) = 
[u ± m]g, cj+j;_j_i = = /^(m) with ft{m) = {[s =F m]q[s ±m + l]q)K yielding a new 

series of integrable models, namely q-spin (4s + 2)-vertex models for different values of s (FIG. lb)). 

2. Another constant but different set of choices: = = 1, Cg" = = giving = ztv^itl 
yield an exponentially deformed Lie algebra [S'+,5~] = ^^sina The representation (^) is valid now 

with g{u) = - — ^^.^^ ' and reproduces the lattice version of the quantum Liouville model. One 
can as well construct the related vertex model with the BW obtained as particular reductions of ( |3.8| ) 
similar to the above case. 

3. On the other hand another constant choice = C2 = 1, = —iq, C2 = ^ leading to 
= 2sina, = 2cosa reproduces the well known g-boson algebra [^1-^] = ^) [^^)-^] = 

-A-^, [A, At] = ^"'^"if^^^^" from by denoting S+ = pA, S' = pA^ = -iV, p = (cot a) 5 
and corresponds to a new quantum integrable q-ho sonic model with the Lax operator reduced from 
(O) as 

' ^ ^ V e-'^^[\+{N +^)]q y ' ^ ' 

Realizing the g-boson through standard boson one may construct an integrable quantum model, 
representing a lattice version of the quantum derivative nonlinear Schrodinger equation (QDNLS) and 
at the continuum limit the corresponding field model. Fusing two such models one can build further 
a quantum integrable massive Thirring model Q. The QDNLS is also related to the exactly solvable 
interacting bose gas model with derivative 5-function potential. 

For constructing the related vertex model we require matrix representation of the q-bosonic opera- 
tors, which with the present choice of the central elements and assuming k = s = 0, n = — m, may be 
derived from (|3?7| ) as < n\N\n >= n5n,n, < n\A'^\n >= fQ{n)6n+i,n, < n\A\n >= ffl{n)5n-i,n 
with f^{n) = ([1 + nUi-n - 1]],)^ = ^jll + ^^ija, foi^) = /-(«-!) = ^[n]J. Consequently the 
BW of the related statistical model reduced from ( ^ ) take the form as 

^±j;±j(n) = ie±^"*[n T (j + - 1)]„ = = /+(i - 1) = -^[j - (4.11) 

with (j) = ^(1 + ^) which generates new q-boson (4n + 2)-vertex model (FIG.lc)). 

4. It is easy to see that for all the following parameter choices: i) = 1 , o = 1, 2, or it) cf = 
±1, or iii) = 1, with rest of the c's being zero, we get = with the underlying algebra 
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reducing to [5+, 5"^] = 0, [S''^,5^] = ±5'='= and ( p.Sp simply to = —ip, = ae^" . This yields 
the relativistic quantum Toda chain recovering also its different Lax operator constructions. 

The statistical systems related to this case however seem to be give uninteresting vertex models 
with infinite matrix traces. 



4.2 Models with q roots of unity 

It is obvious that the integrability of a system remains intact even if q is chosen as solutions of q^ = 1. 
with parameter a taking discrete values Oa = 27r^, a = 1,2, . . . ,p — 1 |jl2|. However this opens up an 
excellent possibility for regulating the dimension of the representation of the underlying algebra. All 
quantum models constructed above will yield their restricted versions when q is taken as the roots of 
unity. Among these models only a few appear to have been studied earlier, which include restricted 
sine-Gordon model connected with the q-spin and restricted DNLS associated with the q-boson. 

The related vertex models generate now completely new series, since the number of the possible 
configurations can be regulated through integer p. To analyze this fact we focus on the action of S~ 
in (|3.7|) assuming k = : and observe that due to = sinOaP = 0, unlike generic q we can get now 
S~ \ — s,s >= at s = p — (s + 1), which reduces matrix (3/7) to a finite dimensional representation. 
As a consequence, for q-spin with fixed p,0 < p < 2s + 1, we get now p — 1 number of different 
{4p — 2)-vertex models for different discrete values of a^. 

The situation becomes more interesting when applied to the q-boson with finite p, since now in place 
of its unbounded representation we obtain a finite dimensional matrix representation leading to an 



intriguing series of (2p— 2)-vertex models with BW described by the same form (4.11) as for the generic 
q-boson case, but now with different possible discrete parameter values q = e*"", a = 1, 2, . . . ,p — 1. 



4.3 Rational class 

At the limit a — > or g ^ 1 the Lax operator ( |3.3[ ) , the BM ( |3.8[ ) as well as the i?-matrix reduce to 
their corresponding rational limits with the underlying algebra becoming 

[s+,s~] = 2m+s3 + m", [s^ s=^] = ±s=^, [m^,-]=0, (4.12) 

with the central elements — > m^. We can find as before various realization of this algebra directly 
from ( |3.5D using however a limiting form of (|3.6D as go{u) = i{{s — u){m^{u + s + 1) + m~))2 and 
construct in the similar way as above different integrable quantum and statistical models belonging 
to this rational class. 

1. Thus at m+ = l,rn-~ = 0, when one gets the standard su{2) algebra, the quantum models 
produced are the XXX spin chain and the lattice NLS model, while the corresponding integrable 
vertex models related to the spin-s operators recover those obtained earlier through fusion technique 
i- 

2. A complementary choice m"*" = 0, m~ = 1, on the other hand yields the bosonic algebra generates 
another quantum simple lattice NLS model |l^ . The corresponding bosonic vertex model represents 
a nontrivial integrable statistical model, apparently never studied before. 
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3. A trivial choice = reproduces the algebra related to the relativistic Toda chain consid- 
ered above, which however in the rational limit gives its well known nonrelativistic variant. The 
corresponding vertex model as before seem to be not physically interesting. 



5 Inhomogeneous and hybrid models 

An immediate generalization of all the above models is possible by considering the central elements 
c's appearing in the Lax operator ( |3.3D to be site dependent functions. This would lead to a new class 
of integrable inhomogeneous extensions of the above models and may be interpreted as models with 
impurities, varying external fields, incommensuration etc. 

Thus we can construct for example variable mass sine-Gordon, Liouville models, variable coefficient 
NLS , Toda chain models etc similarly the corresponding integrable q-spin and q-boson vertex models 
in varying external field. 

Another intriguing class of models can be formed by regulating the inhomogeneity such that the 
Lax operators of different models are arranged to sit at different lattice sites along the chain defining the 
transfer matrix of the model as r(u) = tr (jJpYlf^l^ L^p {u)^ , where (u) indicates Lax operators 
of different models belonging to the same class with the same i?-matrix and occurring times in 
the total number of sites N = J2f3 

Thus one can construct new series of exotic integrable models like hybrid sine-Gordon-Liouville 
model, hybrid NLS-Toda chain or spin-boson model, describing different types of nonlinear interactions 
at different domains of the coordinate space. Similarly one obtains hybrid integrable statistical models 
like (un-)deformed spin-boson models etc. combining different types of vertex models belonging to 
the same class (FIG.l). 

Considering the continuum limit A — > we recover the quantum field models from their respective 
lattice versions constructed above, while their field Lax operator £(x,A) takes the form Lj{\) — > 
/ + iA£(x, A) + 0(A^). The associated i2-matrix however remains the same as its discrete counterpart. 
At the limit /i ^ one obtains the corresponding classical models with the i?-matrix reducing also 
to its classical form R{X) = I + hr{X) + 0(T)?'). For statistical models on the other hand the most 
relevant is the thermodynamic limit: oo, with A, Ti fixed. 



6 Unified solutions 

: There is a well formulated algebraic Bethe ansatz method for exactly solving the eigenvalue problem 
of the transfer matrix t{u) = tr(J^ Li(u)), when the L, R matrices are given [11]. Therefore for our 



unified model represented by the generalized Lax operator ( |3.3D and the i?-matrix ( |2.2| ) we get 

n n r _|_ 11 

A(A) = (< 0|Lii(n)|0 >)^(n) n f{uk -u) + {< 0|L22(u)|0 >)^(n) [] f{u - u^), f{u) = 



k k L-jy 

(6.13) 
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with all possible solutions of {u^} to be determined from the Bethe equations 

/«0|LnMO»\" + ^ , 

Here |0 > is the pseudo- vacuum and the only model dependent parts in both the above equations are 
given by the actions of the upper and lower diagonal operator elements L^^, L^^ of the Lax operator. 
Note that since in our scheme the i2-matrix is the same for all models and their L-operators are 
given through various reductions of the eigenvalue form ( |6.13 ) together with the Bethe equation 
( |6.14D give a unifying scheme for exactly solving all the integrable quantum as well as statistical models 
constructed here. 

For our vertex models the Lax operator elements in the above equations should be replaced by their 
matrix representations expressed through the BW (|3^ ) as (< 0|L^^(u)|0 >) = ix'+_i;+_i(< 0|L^^(u)|0 > 
) = The total number of independent solutions for the eigenstates should be equal to 

the dimension of the vector space on which the transfer matrix acts. For the vertex models in our 
construction this is -ftT = and the partition function should be given through exact eigenvalue 
solutions as Z = limM,N^ooti"v{T^'^) = limM,N^ooY^k ^¥ i ™ the thermodynamic limit. Note that 
though in general the dimension of the transfer matrix may even be infinite (including degeneracies), 
the corresponding partition function must be well defined. The Hamiltonian of the quantum models 
related to such statistical systems would generally involve nonlocal interactions which however are not 
relevant for the associated vertex models we are concerned with. Though the Bethe equation gives the 
form for deriving exact solutions, only in the thermodynamic limit one usually expects to find such 
solutions by converting this algebraic equation into an integral equation For our unified model 
assuming = ^(c^)*, for physical reasons we can similarly derive from ( |6.14| ) 

V{ctu) = 2npiu) - I dvp{v) -. (6.15) 

J cos 2as — cosh{u — v) 

It is important to note that the rhs of the above equation totally coincides with that for the 6- 
vertex model while the Ihs represents the model-dependent part and is expressed through BM 
= c5^e"+*'^" + c5-e-("+^^") = riuj{u) given by (U) as V{cf,u) = '^''^'\~y' ■ Therefore the 
integral equations related to all vertex models including the q-spin and q-boson vertex models can be 
obtained from the single equation ( |6.15 ) at proper choices of cf, which should also naturally include 
the well known case of 6- vertex model Hybrid models, models with q roots of unity or at g ^ 1 
can also be covered by the above unifying scheme or its extensions. 



7 Concluding remarks 

Thus we have prescribed an unifying scheme for constructing as well as solving integrable quantum 
and vertex models of certain classes , which covers lattice and field models of (non-)relativistic and 
(isotropic) anisotropic types as well as the corresponding integrable statistical models. Along with the 
known models one can construct new models, including a novel series of q-spin and q-boson vertex 
models related to q at roots of unity. Inhomogeneous and hybrid models constitute new exotic classes 
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of integrable models. Using algebraic Bethe ansatz we can also find exact solutions for all these models 
in a unifying way. For some other details concerning our unified scheme the readers may consult the 
related works |l5|, 0]. 
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FIG. 1. Integrable vertex models with horizontal (h) links taking 2 values, while the vertical (v) ones may 
have D possible values, a) 6-vertex b) q-spin vertex and c) q-boson vertex models. Combining a,b,c) an 
integrable hybrid model may be formed. qP ^ I gives Z? = p in b) and c) 
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